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LONG-TERM  GOALS 

In  the  context  of  EuroSTRATAFORM,  the  long-term  goal  of  this  and  a  closely  related  project  (“The 
Role  of  Negative  Buoyancy  in  the  Morphodynamics  of  River  Dominated  Shelves:  A  Global 
Synthesis”)  is  to  contribute  toward  simultaneously  understanding  (i)  short  term  oceanic  processes  that 
erode,  transport  and  deposit  sediment  in  the  margin  system  and  (ii)  the  creation  of  the  preserved 
stratigraphic  architecture,  seafloor  morphology  and  sediment  facies  on  continental  margins.  In  order  for 
models  at  such  disparate  time-scales  to  interact,  they  must  communicate  through  expressions  which 
upscale  the  underlying  physical  processes.  Under  Modeling  Task  D5,  the  EuroSTRATAFORM  white 
paper  specifically  states:  “Coherent  techniques  will  be  developed  for  upscaling  individual 
processes/events  into  long-term  stratigraphic-architecture  and  seascape-evolution  models.” 

Another  conclusion  common  to  recent  ONR  Coastal  Geosciences  programs  is  recognition  of  the 
dominant  role  played  by  episodic,  high  energy  events  in  driving  both  sediment  transport  and  bed 
formation.  During  STRATAFORM,  the  majority  of  across-shelf  sediment  flux  was  found  to  be 
associated  with  a  few  major  flood  and  storm  events  occurring  over  just  one  or  two  weeks  every  few 
years.  In  its  section  on  modeling,  the  EuroSTRATAFORM  white  paper  similarly  states,  “Of  special 
concern  will  be  important  events  of  strata  formation  (e.g.,  debris  flows,  extreme  floods)  that  are 
difficult  to  observe.”  Thus  another  long-term  goal  of  both  this  “upscaling”  project  and  its  allied  “global 
synthesis”  project  is  to  specifically  understand  the  role  of  energetic  sediment  transport  in  depositing 
sediment  on  margins  and  shaping  morphology. 

OBJECTIVES 

Our  work  in  FY2004  combined  a  no-cost  extension  of  this  “upscaling”  project  along  with  the  first  six 
months  of  a  closely  related  project  on  “The  Role  of  Negative  Buoyancy  in  the  Morphodynamics  of 
River  Dominated  Shelves:  A  Global  Synthesis”.  The  common  objectives  of  these  two  projects  include 
(1)  investigating  the  role  of  high  energy  sediment  gravity  flows  in  shelf  deposition  at  multiple  time- 
scales,  and  (2)  investigating  the  role  of  these  flows  in  shelf  deposition  at  diverse  spatial  scales  and  in 
diverse  settings  world-wide. 

To  help  address  these  objectives  at  the  shortest  time-  and  space-scales,  we  are  applying  a  wave¬ 
resolving,  high  vertical  resolution  numerical  model  with  advanced  turbulence  closure  to  simulate 
energetic  sediment  suspensions.  The  consistency  of  the  internal  structure  revealed  by  these 
simulations  and  the  bulk  formulations  applied  in  our  wave-averaged,  depth-integrated  analytical 
approach  is  being  examined.  Scales  are  also  being  bridged  from  wave-averaged  (or  “burst-averaged”) 
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quasi-steady  conditions  at  a  single  location  to  deposition  through  whole  flood  seasons  distributed  over 
entire  subaqueous  deltas.  An  important  question  is  whether  bulk  formulations  which  work  over  a 
period  of  hours  also  work  over  periods  of  weeks  to  months.  Finally,  analytical  and  simple  numerical 
approaches  are  being  used  to  upscale  these  same  relations  for  scenarios  where  accretion  feeds  back  to 
impact  hydrodynamics.  The  largest  scales  compare  diverse,  globally  distributed  environments  and 
incorporate  subsidence,  tectonism  and  changes  in  sea  level. 

APPROACH 

Seconds  to  hours:  To  provide  guidance  in  upscaling  realistic  turbulence  closure  models  under  energetic 
sediment  transport  conditions,  we  are  running  numerical  model  simulations  of  sediment-stratified 
bottom  boundary  layers  using  the  1-D  General  Ocean  Turbulence  Model  (GOTM,  Burchard,  2002). 
GOTM  is  an  open  source,  FORTRAN-based,  two-equation  (k-epsilon)  turbulence  closure  model 
which  uses  the  most  recent  formulations  for  stability  functions  shown  to  perform  well  in  the  presence 
of  strong  thermohaline  stratification.  We  have  added  suspended  sediment  to  this  model  and  are 
investigating  model  behavior  as  one  approaches  the  limit  of  critical  sediment-induced  stratification 
(Scully  and  Friedrichs,  2003).  We  are  now  systematically  examining  the  parameter  space  of  varying 
sediment  settling  velocity  (ws)  from  hindered  settling  in  fluid  mud  (ws  <  0. 1  mm/s)  to  sand  (ws  >  10 
mm/s)  and  the  wave  period  of  the  bottom  boundary  layer  from  storm  waves  (order  10  s)  to  tides  (order 
10  hours).  We  are  comparing  these  simulations  to  available  field  and  laboratory  data  of  very  high 
concentration  fluid  mud  layers  within  momentum  deficit  layers  under  waves  (e.g.,  Eel  shelf,  Traykovski 
et  al.,  2000)  and  tides  (e.g.,  Amazon  shelf,  Trowbridge  and  Kineke,  1994)  and  of  critically  stratified 
constant  stress  layers  of  sand  under  waves  (e.g.,  Hanover  Wave  Flume,  Dohmen-Janssen  and  Hanes, 
2002)  and  of  flocculated  mud  under  tides  (e.g.,  York  River,  Friedrichs  et  al.,  2000). 

Hours  to  months:  The  analytical  relations  governing  a  wave-supported  sediment  gravity  flow  in  a 
wave-averaged,  vertically-integrated  sense  are  (Wright  et  al.,  2001;  Scully  et  al.,  2002): 

a  g  s  C  ps  1  =  cd  lul  ug,  Ri  =  g  s  C  ps_1  lul  2  (1,2) 

In  Eqs.  (1)  and  (2),  a  is  the  sine  of  the  bed  slope,  g  is  the  acceleration  of  gravity,  ps  is  the  density  of 
siliceous  sediment  and  s  is  its  submerged  weight  relative  to  sea  water,  C  is  the  depth-integrated  mass 
concentration  of  suspended  sediment  within  the  wave  boundary  layer,  ct|  is  the  bottom  drag  coefficient, 
and  Ri  is  the  gradient  Richardson  number.  The  key  velocities  associated  with  Eqs.  (1,2)  are  the 
downslope  velocity  of  the  gravity  current  (ug),  the  rms  amplitude  of  wave  orbital  velocity  (uw)  and  the 
absolute  amplitude  of  the  instantaneous  velocity,  |u|  ~  (uw“  +  ug")  ",  all  evaluated  near  the  top  of  the 
wave  boundary  layer.  By  combining  (1)  and  (2),  one  can  solve  explicitly  for  ug  and  C  and,  thus,  for  the 
downslope  sediment  transport  rate,  Qg  =  ugC,  and  the  deposition  rate,  D  =  -  dQg/dx.  It  is  relatively 
straightforward  to  force  these  equations  with  observed  wave  heights  and  river  discharge  to  predict 
gravity  flow  speed,  concentration  and  deposition  over  observed  bathymetry.  This  has  been  our 
approach  in  modeling  gravity  flows  over  hours  to  entire  flood  seasons  for  STRATAFORM  and 
EuroSTRATAFORM  field  sites  (Wright  et  al.,  2001;  Scully  et  al.,  2002,  2003;  Scully  and  Friedrichs, 
2004;  Friedrichs  and  Scully,  2004b,  in  prep.). 

Years  to  centuries:  We  are  taking  two  approaches  in  upscaling  from  a  season  to  decades  and  beyond. 
The  first  approach  uses  predicted  deposition  (plus  any  appropriate  shifts  in  base  level  or  sea  level)  to 
change  the  bathymetry  and,  thus,  the  hydrodynamics.  This  quasi-numerical,  long  time-scale  approach 


follows  the  same  basic  philosophy  as  the  hours-to-months  simulation,  but  using  annual  representative 
wave  heights  and  discharge,  and  with  bathymetric  feedback  affecting  subsequent  hydrodynamics 
(Friedrichs  and  Scully,  2004a, b,  in  prep).  Our  other  long  time-scale  approach  looks  for  equilibrium 
analytical  solutions  to  bathymetric  profiles  consistent  with  the  physics  of  wave-supported  gravity 
flows.  One  option  in  this  regard  is  to  define  an  equilibrium  shelf  profile  as  that  which  causes  gravity 
flows  to  be  just  strong  enough  to  carry  the  riverine  sediment  supply  into  deep  water  without  any  net 
deposition  (Friedrichs  and  Wright,  2004).  Conceptually,  this  applies  well  to  a  situation  where  a  delta 
has  prograded  and  steepened  to  the  point  that  the  wave-supported  gravity  flows  now  just  bypass  the 
shoreward,  concave  downward  portion  of  the  delta,  producing  no  further  deposition  in  this  region. 
Another  option  is  to  define  as  equilibrium  a  sigmoidal  clinoform  or  subaqueous  delta  with  a  fixed 
profile  which  progrades  across  a  flat  paleo-shelf  at  a  constant  rate  (Friedrichs  and  Scully,  2004a, b,  in 
prep).  These  analytical  solutions  for  shelf  profiles  employ  additional  continuity-related  equilibrium 
conditions  beyond  Eqs.  (l)-(2). 

WORK  COMPLETED 

In  FY2004  we  continued  the  process  of  publishing  our  related  work  from  the  ONR  STRATAFORM 
project  (Hill  et  ah,  2004;  Parsons  et  ah,  2004).  Our  newest  research  efforts  in  FY2004  focused  mainly 
on  the  time  scales  of  hours  to  centuries.  Over  scales  of  hours  to  months,  we  applied  Eqs.  (l)-(2)  to 
seasonal  flood  deposition  in  the  region  of  the  Po  subaqueous  delta  (Figure  1)  (Scully  and  Friedrichs, 
2004;  Friedrichs  and  Scully,  2004b,  in  prep.).  Over  scales  of  years  to  a  century,  we  applied  the  same 
equations  to  the  evolution  the  Po  delta  since  the  mid- 1 800 ’s  acceleration  in  progradation  of  the  Pila 
distributary  (Figure  2)  (Friedrichs  and  Scully,  2004a, b,  in  prep.).  In  the  context  of  long-term 
equilibrium  shelf  profiles,  we  finished  our  first  comparison  of  equilibrium  shelf  profiles  around  the 
world  associated  with  bypassing,  wave-supported  gravity  flows  (Friedrichs  and  Wright,  2004;  Pratson 
et  al.,  2004a, b).  We  also  derived  analytical  solutions  for  similarly  force  equilibrium  prograding 
clinoforms  and  applied  the  solution  to  the  Po  shelf  (Friedrichs  and  Scully,  2004a, b,  in  prep.). 

RESULTS 

Figure  1  displays  forcing  conditions  and  flood  deposition  around  the  Po  subaqueous  delta  as  predicted 
by  Eqs.  (l)-(2)  when  applied  to  the  large  flood  of  fall  2000.  For  the  region  encompassing  the  Po 
subaqueous  delta,  the  equations  governing  wave-supported  gravity  flows  were  implemented  in  the 
same  manner  as  that  described  for  the  Eel  shelf  by  Scully  et  al.  (2003).  Bottom  orbital  velocity  was 
determined  by  linear  wave  theory  applied  to  wave  heights  observed  offshore  of  Ancona,  Italy. 

Sediment  discharge  from  the  Po  River  before  separation  into  its  distributaries  was  initially  estimated  by 
applying  the  rating  curve  of  Syvitski  et  al.  (2004)  for  a  single  river  channel  neglecting  random  variability 
(Figure  la).  Sediment  input  was  reduced  to  56%  of  that  predicted  by  Syvitski  et  al.  (2004)  and  further 
restricted  to  the  second  half  of  the  flood  when  discharge  and  waves  occurred  together.  The  along-shelf 
distribution  of  sediment  used  to  represent  the  various  distributaries  is  indicated  by  the  dashed  red  line 
in  Figure  lb.  Results  are  consistent  with  the  general  distribution  of  the  flood  layer  as  mapped  by 
EuroSTRATAFORM  investigators  and  suggest  wave-supported  gravity  flows  can  account  for  the 
majority  of  the  observed  deposit.  In  particular,  modeled  and  observed  deposition  both  occur  in  deeper 
water  off  the  Pila  than  off  the  Tolle,  and,  in  both  cases,  maximum  deposition  occurs  near  the  region  of 
maximum  bed  slope. 

Figure  2  displays  results  of  a  long-term  simulation  of  subaqueous  delta  evolution  using  Eqs.  (l)-(2). 


The  simulations  mimicked  profile  evolution  off  of  the  Pila  and  Tolle  distributary  mouths  since  the  mid 
1800’s.  Each  case  started  with  the  same  linearly  sloping  initial  bathymetry  and  was  subjected  to  the 
same  1  m  wave  height  and  the  same  rate  and  spatial  distribution  of  subsidence.  The  only  difference 
was  the  rate  of  sediment  input,  which  was  five  times  larger  for  the  Pila  simulation  than  for  the  Tolla 
case.  For  the  Pila  case,  deposition  overwhelmed  subsidence,  and  the  subaqueous  delta  prograded 
seaward  (Figure  2a),  reaching  a  shape  and  progradation  rate  very  similar  to  our  analytical  equilibrium 
solution.  For  the  Tolle,  in  contrast,  subsidence  overwhelmed  deposition,  and  the  shoreward  portion 
delta  receded  landward  and  steepened  (Figure  2b).  These  results  are  largely  consistent  with  observed 
evolution  of  the  Po  subaqueous  delta  since  the  mid-1800’s  as  reported  by  Nelson  (1970).  In  the  final 
year  of  the  simulation,  maximum  deposition  occurred  in  notably  deeper  water  off  the  Pila  than  off  the 
Tolle,  similar  to  the  patterns  seen  in  observations  and  simulations  of  the  fall  2000  flood. 

In  summary,  our  major  results  in  FY2004  applied  our  simple  model  for  wave-supported  gravity  flows 
to  sediment  deposition  at  time  scales  ranging  from  a  single  major  flood  to  steady-state  progradation  of 
an  equilibrium  clinoform.  The  modeling  approach,  which  limits  sediment  load  via  a  critical  Richardson 
number,  is  applicable  to  fine  sediment  transport  near  river  mouths  wherever  wave  energy  is  available  to 
move  abundant  sediment  offshore  during  floods.  Results  suggest  this  phemomena  can  account  for  the 
majority  of  the  fall  2000  flood  deposit  mapped  by  EuroSTRATAFORM  investigators  in  the  vicinity 
of  the  Po  subaqueous  delta  and  also  for  the  rate  of  delta  progradation  observed  off  the  dominant  Pila 
outlet  of  the  Po  over  the  last  150  years.  Model  results  predict  that  convergence  of  down-slope 
sediment  transport  by  wave-supported  gravity  flows  increases  with  bed  slope  but  decreases  with  slope 
gradient,  such  that  greatest  deposition  occurs  near  where  steep  slopes  first  stop  increasing  with 
distance  offshore.  Thus  on  profiles  which  reach  maximum  steepness  near  shore,  like  those  off  Tolle- 
Gnocca-Goro  mouths  today  or  off  the  Pila  mouth  150  years  ago,  gravity-driven  deposition  occurs  in 
shallower  water.  Over  time,  if  deposition  overwhelms  subsidence,  the  subaqueous  delta  becomes  less 
steep  near  shore  and  steeper  offshore,  and  the  locus  of  deposition  moves  progressively  into  deeper 
water.  If  the  subaqueous  delta  is  prograding  across  a  relatively  flat  shelf,  the  form  of  the  prodelta 
eventually  reaches  an  approximate  equilibrium  which  progrades  seaward  as  a  unit. 

IMPACT/APPLICATIONS 

A  present  limitation  in  long-term  modeling  of  continental  margin  evolution  is  realistic  inclusion  of 
hydrodynamic  processes  driving  shelf  deposition.  Based  on  field  observations  collected  over  the  last 
20  years,  complex  wave-averaged  currents  driven  by  winds  and  pressure  gradients  have  been  thought  to 
be  mainly  responsible  for  cross-shelf  sediment  transport  and  flux  convergence  on  energetic  accretionary 
shelves.  Unfortunately,  it  may  be  exceedingly  difficult  to  predict  wind-  and  pressure-driven  near-bed 
currents  with  sufficient  accuracy  to  produce  realistic  deposits  over  geological  time-scales.  The  ONR 
STRATAFORM  and  EuroSTRATAFORM  projects,  however,  have  identified  a  distinctly  different 
mechanism  for  across-shelf  mud  transport  associated  with  gravity-driven  flows  of  fluid  mud  within  the 
wave  boundary  layer.  Gravity  flows  within  the  WBF  can  be  realistically  modeled  based  on  knowledge 
of  fine  sediment  supply,  approximate  wave  height  and  bathymetry  if  one  assumes  that  the  critical 
Richardson  number  within  the  WBF  determines  the  maximum  capacity  of  the  gravity  flow  to  transport 
mud.  Complex,  externally  forced  mean  currents  do  not  appear  to  play  a  critical  role  in  this  newly 
identified  transport  mechanism.  Thus  the  analytical  model  presented  here  has  the  potential  to  greatly 
reduce  the  complexity  and  computational  requirements  presently  limiting  our  ability  to  perform 
realistic  long-term  simulations  of  the  geologic  evolution  of  many  continental  margin  environments. 
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Figure  1.  (a)  Forcing  conditions  for  the  fall  2000  Po  flood  simulation,  (b)  Predicted  thickness  of 
the  resulting  flood  layer.  Predicted  deposition  is  consistent  with  the  general  distribution  of  the 
flood  layer  as  mapped  by  EuroSTRATAFORM  investigators. 


Figure  2.  Evolving  bathymetry  for  simulation  of  subaqueous  delta  deposition  off  of  the  (a)  Pila  and 
(b)  Tolle  distributaries,  (c)  Representative  profile  locations,  (d)  Predicted  deposition  during  the 
final  year  of  simulation  along  with  the  degree  of  annual  subsidence.  Both  the  long  term  and  final 
year  simulations  are  largely  consistent  with  observations  from  the  Po  subaqueous  delta. 


TRANSITIONS 


Our  data  on  bed  stresses  and  resulting  sediment  resuspension  from  earlier  years  of  our 
STRATAFORM  project  have  been  made  available  to  others  and  are  being  used  to  verify  various 
bottom  boundary  layer  and  sediment  transport  models.  Our  data  can  easily  be  accessed  via  data 
reports  (which  include  data  summaries  on  diskettes)  and  via  the  VIMS  STRATAFORM  website: 
http://www.vims.edu/physical/projects/CHSD/projects/ONR/index.html.  Published  papers  by  others 
which  have  directly  utilized  VIMS  data  include  Morehead  and  Syvitski  (1999),  Ogston  et  al.  (1999, 
2000),  Reed  et  al.  (1999)  and  Zhang  et  al.  (1999).  Additional  papers  by  non- VIMS  authors 
incorporating  VIMS  data  are  in  preparation.  Our  analytical  formulation  for  sediment  flux  and 
deposition  by  critically-stratified,  gravity  flows  has  already  been  incorporated  into  long-term 
simulations  of  margin  stratigraphic  development  by  James  Syvitski's  group  (Syvitski  et  al.,  2001, 
2002).  Our  analytical  approach  has  also  been  made  available  to  other  modelers,  such  as  Fan,  Harris, 
Niederoda,  Reed,  Swift,  and  Traykovski,  all  of  whom  are  at  various  stages  of  incorporating  gravity 
flows  into  more  complex  numerical  simulations  of  shelf  sedimentation. 

RELATED  PROJECTS 

The  following  active  projects  involving  Friedrichs  also  focus  on  coastal  sediment  transport: 

1.  Forecasting  Scour  Related  Mine  Burial  Using  a  Parameterized  Model.  Office  of  Naval  Research 
(www.vims.edu/physical/projects/CHSD/projects/MBP). 

2.  Sediment  Dynamics  of  a  Microtidal  Partially-Mixed  Estuary.  National  Science  Foundation 
(www.vims.edu/physical/projects/CHSD/projects/CAREER). 
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